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<T> ! Abstract 

cn : 

We look for the global in time solution of the Cauchy problem corre- 
sponding to the asymptotically flat spherically symmetric E.V.M system 
with small initial data. Using an estimate, we also prove that if solution 
of the system stated above develops a singularity at all time, then the first 
one has to appear at the center of symmetry. 
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The EVM system models the time evolution of self-gravitating collisionlcss 
charged particles in the context of general relativity. The particle could be 
for instance the electrons in a plasma. As it is proved in [Hj , that system in the 
^jy context of spherical symmetry has nice properties. Firstly, the electromagnetic 

field F a /3 created by the fast moving particles reduces to its electric part and 
secondly if / = then F a p = and the latter is not true in the case without 
. spherical symmetry, where it is possible to have non trivial source-free solutions 

^ ' of the Maxwell equations. 

In [7] , the authors established with the help of an iterative scheme one local 
existence theorem and a result on the continuation criterion of solutions is de- 
duced. Note that these results are valid only for the Cauchy problem with small 
initial data. Here, we aim to establish a global existence theorem of solutions for 
the asymptotically flat spherically symmetric EVM system. To do so, contrary 
to the uncharged case, we defined a set of initial data in such a way that for 

o o 

fixed /, the solution (A, e) of the constraint equations(see satisfies the con- 

o 

dition that A is bounded in the L°°-norm. The latter is possible, provide that 

° o 

the charge q is small. Note that in the above, A (resp. e, /) is the initial datum 
of the metric function A(£, r)(resp. electric field e(t, r), distribution function of 
the particles /). So, the global existence theorem we obtain, generalize that 
established by G. Rein and AD. Rcndall in |12|. 
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Next, for the related literature, global existence for the Vlasov-Maxwell sys- 
tem without symmetry is not known in general. The most general result which 
has been proved is for a case with one-dimensional symmetry group |3] . Global 
existence for the Yang-Mills equations in Minkowski space was proved by Eard- 
ley and Moncrief [5J and another proof by quite different methods was given by 
Klainerman and Machedon [S]. Global existence for the Yang-Mills equations 
on a general globally hyperbolic spacetime was proved by Chrusciel and Shatah 

m 

Now, in the case of spherical symmetry, the initial boundary value problem 
for the asymptotically flat EVM system admits the global solution and the 
obtained spacetime is geodesically complete i.e each geodesic is defined on R. So, 
it is well known in general relativity that this kind of solution does not develop 
a singularity. Concerning the generic initial data, it is proved by A.D. Rendall 
[T^] that a singularity may occur. Also, there is a related result of G. Rein , 
A.D Rendall and J. A Schaeffer JI] which states that for the asymptotically flat 
spherically symmetric Einstein- Vlasov system, if solutions develop a singularity 
at all time then the first one has to occur at the center of symmetry, and also 
if particles remain away from the center, the local solution of Cauchy problem 
for that system with a given initial data can be extended to obtain the global 
one. That results concern uncharged particles. 

It would be interesting to see what happens in the case of charged particles. 
More precisely, is it possible to generalize the above results in the context of 
charged particles? The answer is yes and this is the second important result of 
this paper, the first one being the global existence theorem we stated earlier. 

The paper proceeds as follows. In Sect. 2, we formulate the Cauchy problem. 
In Sect. 3, we show that if the charge is small then one can choose a set of initial 
data so that solutions depend continuously on the latter. In Sect. 4, we establish 
a global existence theorem of solution for the corresponding initial boundary 
value problem and the obtained spacetime is geodesically complete. In Sect. 5, 
we establish the main estimates we use to prove the global existence theorem 
in the case that all particles remain away from the center. In Sect. 6, we prove 
a local existence theorem and a continuation criterion for the exterior problem. 
Sect. 7 contains the essential result obtained using all what have been established 
in the last two sections above. 

2 The Cauchy Problem 

We consider fast moving collisionless particles with charge q. The basic space- 
time is (R 4 , g), with g a Lorentzian metric with signature (— , +, +, +). In what 
follows, we assume that Greek indices run from to 3 and Latin indices from 1 
to 3, unless otherwise specified. We also adopt the Einstein summation conven- 
tion, both the gravitational constant and the speed of light are setting to unity. 
Now, Using the assumption of spherical symmetry, one obtains with the help of 
the Killing equation for instance that g can be taken of the following form [Sj: 

ds 2 = -e 2 »dt 2 + e 2X dr 2 + r 2 {d9 2 + (sm0) 2 d<p 2 ) (2.1) 



2 



where [i = fi(t,r); A = A(i, r); t e R; r e [0,+oo[; e [0,7r]; <y9 e [0, 27r]. One 
shows (see 8 ) that the spherically symmetric EVM system can be written as 
the following P.D.E system in A, fx, f, e: 



-2A 



e 



(2rA' - 1) + 1 = 8irr 2 p (2.2) 



e _2A (2r/i' - 1) + 1 = 8nr 2 p (2.3) 



dt y/l + v 2 dx 



\ r J r ov 

(2.4) 

( r Ve)=gr 2 e A M (2.5) 

or 

where A' = ^; A = ^ and: 

p(t,x)= / /(t, i, v) y/l + v 2 dv + ~ e 2X ^e 2 (t, x) (2.6) 

p(t,x)=[ (™\\{t,x,v)^^-\e^e\t,x) (2.7) 
Jr3 V r J \/l + v 2 2 

M(t,S) = / f(t,x,v)dv. (2.8) 

Here (12.21 and l|2.3[) are the Einstein equations for A and /Lt, (|2.4II is the Vlasov 
equation for /, / being the distribution function of the charged particles which is 
defined on the mass-shell and l|2. 5fl is the Maxwell equation for e. Here x := (x l ) 
and v belong to I 3 , r :=| x |, x.v denotes the usual scalar product of vectors 
in M 3 , and v 2 := v.v. The distribution function / is assumed to be invariant 
under simultaneous rotations of x and v, hence p, p and M can be regarded 
as functions of t and r. It is assumed that f(t) has compact support for each 
fixed t. We are interested in asymptotically flat spacetime with a regular center, 
which leads to the boundary conditions that: 

lim Aft, r) = lim fx{t, r) = \(t, 0) = lim e(t, r) = e(t, 0) = (2.9) 

r— >oo r—*oo r— >co 

Now, define the initial data by: 

{o o o 

/(0, x, v) = f(x, v); A(0, x) = X(x) = A(r) ^ 1Q . 

fi(0, x) = fi(x) = £(r); e(0, S) = e(£) = e(r) 

o 

with / being a C°° function with compact support, which is nonnegative and 
spherically symmetric, i.e 

VA e SO(3), V(x, v) e K 6 , f(Ax, Av) = f(x, v). 
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We aim to solve the initial boundary value problem H2.2[l . (|2.3II . (|2.4|) . 12. 5|) . 
US and (|2~TT)|) . 

Remark 2.1 7Vo£e £/iat i/ie assumption e(t,0) — makes sense. To see the 
latter, we use the fact that the electric field E(t,x) :— e(t, ^Oifr *s spherically 
symmetric in the x variable i.e AE = E for every A £ SO(3) and once 
E(t, 0) = 0, one easily deduces the announced result. 

3 Continuous dependence on the initial data 

In this section, we prove that solutions depend continuously on their initial data. 
Besides being of interest in itself for physically viable theory ( cf.[|16j. p. 243f]), 
the results of this section will be applied in the proof of global existence for small 
data in the next section. Before doing so, we recall the local existence theorem 
we proved in 7 on which our present result relies. The constraint equations are 
obtained setting t = in 1)2.2(1 . 1|2.3|) and l|2. 5JI . Using perturbation techniques, 

o 

it is proved in 6 that for / small the constraint equations has a global smooth 
(i.e C°°) solution. First of all, we need to make precise the notion of regularity 
we use in the sequel. 

Definition 3.1 A solution (f,X,/j,,e) of Wty . fO|) . and fO)) . is said to 

be regular if: 

(i) f is nonnegative, spherically symmetric, C , and f(t) is compactly sup- 
ported for all t, 

(ii) A > 7 /i < 0, and X, /i, A', // and e are C 1 which X, /i and e satisfying 
the boundary condition \2.9ji . 

o 

Theorem 3.1 (local existence, continuation criterion) Let f £ C°°(M 6 ) 
be nonnegative, compactly supported and spherically symmetric such that: 

8lT Jo (/ -f( r > v ^ 1 + v2dv ) ds<r t 3 ' 1 ) 

Let A,/i,e € C°°(R 3 ) be a regular solution of the constraint equations. Then 
there exists a unique regular solution (A, \x, /, e) of the asymptotically flat spher- 

o ° o 

ically symmetric EVM system with initial data (A, /i, /, e) on a maximal interval 
I CM of existence which contains 0. If 

sup{| v I I (t,x,v) £ supp/, t > 0} < +oo 

then sup / = +oo, if 

sup{| v | | (t,x,v) £ supp/, t < 0} < +oo 

then mi I = — oo. 
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Proof: See [7|. 

For ro > 0, Mo > and A > 0, we consider the following set of initial data: 
D := {(/,A,e) £ C*°°(R 6 ) x (^([O, +oo[)) 2 , / > 0, spherically symmetric, 

° ° o 

and satisfies (I3.1[l supp/ C B(ro) x B(uq) and (A,e) 

o 

is a regular solution of the constraints with || A ||l°°< A}. 

Fix the solution (X g , [i g , e g , g) of (|2.2|) . I|2.3|) . (|2.4[) and (|2.5fl with initial datum 

(X g ,e g ,g) £ D and right maximal existence interval(r.m.e.i) [0, T g [. We aim to 
control the distance of another solution (A/,^/,e/,/) from (A g , /i g , e g , g) and 
the relation between the maximal existence times Tt and T g in terms of distance 
of the initial data, [0, T/[ being the r.m.e.i of (A/, fj,f, e/, /); the whole argument 
would also work for t < 0. To do so in our context, we first have to control the 

distance between two solutions (A/, e/), (X g , e g ) of constraint equations and the 

essential tool we use is the fact that we can construct a set of initial data (A, e) 

o 

such a way as the L°°-norm of A is uniformly bounded. This comes from the 
continuous dependence of solutions of the constraint equations on parameter q, 
when q is small as it is shown in 0. 

3.1 Continuous dependence of solutions for the constraint 
equations 

Let us give the main result of this section: 

o o o ° o 

Proposition 3.1 Consider (f,\f,et),(g,\ g ,e g ) £ D. Given a sufficiently 
small real number e > 0, if d :=|| / — g ||z,°°< £, then 

IIA/IU-, \\e f U~<C (3.2) 



(3.3) 



||e A /-e A «|U~, \\e x fe f -e x 4 g \\ L ~, || 1/ - I fl |U- 
IIA/-AJU-, || e A f°e) ~e 2 ^t g \\ L ^, \\ e A ' - e 2 ^ || L ~< Cd 

o 

where the constant C depends on ro,Uo,3, A and not on f. 

o o o 

Proof : Take (/, A/, e/) £ D, with d < e. The bound of A/ comes immediately 

o 

from the definition of D. Next, setting t — in l|2.5|l and replacing / by /, we 
obtain by integration: 



° ef (r) = ± e -° x f^ f s 2 e° x f^ M f (s)d, 
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where M / is defined as in (|2.8|l . replacing / by /. Distinguishing the cases 

r < tq and r > r$, we obtain, since / is with compact support, | e(r) |< C, 
where C = C(ro,uo,g, A) is a constant , and (|3.2|) holds. To prove inequalities 
(jjOJ), one writes: 



° ° 9 ° ° m r 



o o 

' 1 _l_ c^a— Af 



1 + e A s _A / 

where my = m/(0,r), m(t,r) being the mass function: 

m(t,r) :=4n [ s 2 p{t,s)ds (3.4) 
Jo 

and obtain using the Gronwall lemma: 

| e x f -e A » | (r) < C || f-g || L « exp (c 1 ^ (r - s) 2 | M/(s) | ds^) < Cd 
We also deduce, using the mean value theorem and for e sufficiently small: 
\X f -X g \ (r)=|Loge^-Loge^ | (r) 

< sup jj| \z-e Xs | (r) < CeJ | e A ^ - e A » | (r) 

< Cd, 

and the rest of inequalities in 1|3.3|1 follow immediately from those above. Thus, 
this ends the proof of Proposition 3.1. 

3.2 Continuous dependence of solutions for the Einstein- 
Vlasov-Maxwell on initial data 

We now state the essential result of this section: 

Theorem 3.2 There exists a constant e > 0, a positive increasing function 
£ € C([0, T g [), and a positive decreasing function a <E C(]0, e[) such that 

lim a(f3) = T q (3.5) 

/3->0 

° o 

and for any solution (A/, fif, ef, f) with the initial datum (/,A/,e/) £ D satis- 

o 

fying d :=|| f — g \\l°°< e, we have the estimates: 

Tf > a(d) (3.6) 



II /(*) - 9(t) |U- + II A/(t) - A ff (t) + || /x/(t) - /i 9 (t) + || e f (t) - e g (t) \\ L ~ 
+ || e 2X f^ - e 2A,(t) || ioo + || ^(i) _ A ff (t) 

+ ii n' f (t) - ^(t) \\ L ™<md 

(3-7) 

/or i S [0, o"(d)]. TTie analogous assertion holds for t < 0. 
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Proof : To establish the estimates for the difference to two solutions at time 
t > we first determine a time interval on which we get a uniform bound on 
the supp/(£) for all solutions which the initial data are in D and close enough 
to g. Define 



T (f) : = sup{i €}0,min(T f ,T g )[\ such that < s < t, 

II e f (s)-e g (s) lUoo+H e^e)-e^e\ (s) 
II A/(s) - X g (s) |k- + || /4( S ) - ^(.s) |k-< 1}. 

Notice that for d small enough, say <i < E\ for a suitable defined E\ > 0, the 
estimate defining T (f) holds at t — so that by continuity T (/) > 0. For a 
solution (/, A/, /i/, e/) with e? < E\ the following estimates for the characteristics 
hold on [0,T (/)[: 

|i(«)l<l, 

I *(*) I < A/(«) |U- + 1| ^(s) Hi- + || e/ ( s ) |k~)(i+ 1 «(*) I) 

< C(l+ || A/(s) |k- + || fj,' f (s) Hi- + || e/ (») |k-)(l+ I w(«) I) 
with C = C(g) being a constant. Via the Gronwall inequality this implies that 

supp/(t) c {(x,v) eM 6 | I sc |< r +t, | v \< U g (t)} (3.8) 
for t G [0,T (/)[, where 



+ II /*» Ik- + II e fl («) |k«)ds 



C^:=(l + uo)exp C / (1+ || A ff ( 



Let C be a continuous, increasing function on [0, T g [ which depends only on 
{g, X g , /ig, e g ). Then we obtain the following estimates on [0, Tb(/)[: 

II Pf(t) ~ Pg(t) Hi- + II Pf(t)-Pg(t) |k- + || - fc 9 (t) Ik- 

+ || M f (t) M g (t) || L -< C(f)(|| /(*) - <?(*) |k- + || e^e}(t) - e 2X ^e 2 g (t) |k-) 

(3.9) 

where the matter quantity 



k(t, x) 



x.v 



f{t,x,v)dv 



comes from the following part of the Einstein equations : 

X(t,x) = -Anr(e x+Ii k)(t,x) 

Now, we have: 



(3.10) 



(3.H) 



m f (t,r) m g (t,r) 



m f (t,r) m g (t,r) 



< C(t)(\\ f(t) - g(t) |k- + || e^e){t) ~ e 2X ^ei(t) |k-)- 
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The latter is obtained using the following equality deduced from the integration 
of (EjJoii [r +t,r}: 

e x ^e{t, r) = (j^y^j e x ^' ro+t h(t, r + 1), re [r + t, +oo[. 

Next as the second step, we derive an estimate for the time evolution of / — g 
along the characteristics Zf = (Xf,Vf) corresponding to / (see [7], Proposition 
3.1). Using the fact that / is constant along these characteristics, the mean 
value theorem and since fi — A<l,/x + A<l both for / and <?, one obtains the 
following estimate: 



ds 



(/ - g)(s, Z f (s,t, z)) < C(s)(\\ X f (s) - X g (s) \\ L ~ + || n f (s) ~ fi g (s) \\ L ~) 

+ C(s) || e f (s)-e g (s) \\ L o, 

+ C(s)(\\ X f (s) - X g (s) |U- + II - »' g (s) ||z~); 

and integration of the above w.r.t s on [0, t] yields: 



II /(*) - g{t) ||z- <\\f-9 Hz- +C(t) / || X f (s) - X g (s) || z . ds 



C(t) / || fx f (s) - M9 (s) ||z» <k 



C(t) / (|| e/( S ) - e g (s) + || - A ff (s) || L «)d a 



+ / ii a*>(«) - /*;(«) iu- ds. 

Jo 

(3.12) 

Using l|3.9|l and the fact that || fit) ||z— = || / ||z— <|| 9 ||z- we obtain the 
estimates 



II M*) IU- + || P/(<) Hz- 

+ || k f (t) || L „ + || jv/(*) |U-< c(t), t < Mf) 



(5.10') 



and 



m/(i,r) 



< C(t), t<T f , r>0. 



where the quantity 



N(t,x) :-- 



x.v 



-f(t,x,v)dv 



comes from a following part of the Maxwell equations |7] : 



(3.13) 



(3.14) 
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Using once again the mean value theorem, the inequalities above and (|3.11|) one 
obtains: 

II A/(i) - X g (t) || L » < C(t)(\\ /(*) - </(*) Ik- + II e^e%t) - e^e%t) \\ L ~) 

+ c(t)(\\ \f(t) - x g (t) \\ L oo + || (if(t) - n g {t) Hi-). 

(3.15) 

Now, (|2.3|) yields (—re 2A )' = 8irr 2 p and the integration of this on [0, r] gives: 

e-^) = l- 2m(t ' r) , (3.16) 
r 

So, combining (|2.3I) and (|3.16|l . one obtains: 

fl '(t,r)^e 2 ^\ 1 ^^-+Anrp(t,r)) (3.17) 

and using i|3~T7|l . 1(33)1 and the fact that 

lim /j,f(t,r) = lim / u s (t,r)=0, 

r — >+oo i — > + oo 

we can deduce by integration on [0, +oo[ w.r.t r the following inequality: 

II M/M n a {t) \\ L - < C(t)(|| /(*) - |U- + II ^ (t) - e 2A ^> 
+ C(t) || e 2X ^e 2 f (t)-e 2X ^e 2 g (t) \\ L - . 

(3.18) 

We have, since A + ^ < for both / and g: 

| d s e 2X f - d s e 2X ° | (s,r) < C(s)(\\ e 2A /W _ e 2A s ( s ) || £oo + || A/(s)-A 9 (s) Hz-). 
The insertion of H3.15|l using (|3.18|) in the inequality above gives, by integration: 

|| e 2A /(t ) _ e 2A 8 (t) <|| /X, _ /X g || i<jc +c(i) /"* n _ g{s) dfi 

JO 

+ C(t) / || A/^-A^a) <2s 



Jo 

+ C(t) / || m/(s)-M s (s) IU~)ds 
Jo 

+ C(t) / || e 2A ' (s) -e 2A * (s) Hi- ds 
Jo 

+ C(t) Ae^e^-e 2 ^ 2 ^) |U» ds. 
Jo 

(3.19) 

Now, to obtain an estimate for the last term on the right hand side of (|3.19() . 
we use (|3.14f) and the fact that since /j, < and A + fi < 0, one gets: 

\N g (s,r) \<(r +t) \\g\\ L -<C(a) 

| e M/ _ e M s | ( Sjr ) <| M/ _ Mff | ( S)T -) (Mean value theorem), 
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to deduce by integration that: 

|| e^e)(t) - e^el{t) <|| eM) - e^°e g ||; 



+ / C(s) || f(s)-g(s) Hi- da 



+ / II -M s (s) IIl- rfs 

Jo 

+ fc(.) || e A 'W e/ ( a )_ e *«M efl ( fl ) I), 



(3.20) 

To deduce an estimate for the last term in the right hand side of (|3.20() . we use 
once again l|3.14J) and obtain by integration: 

|| e VW e/ (i) _ e x ^e g (t) || £ - <|| e l 'e f - M g || £ - 

+ ^ C( S )||/( S )- 5 ( S )|U»d 5 (3.21) 

+ / C(«) II M/( s ) - P-g{s) ds. 
Jo 

So, using l|3.21[l . (|3.20|) . i|3.19[l . I|3.12l) and the Gronwall inequality, one obtains, 
since C £ C([0, T 9 [) depends only on 5 and can be taken strictly increasing with 
lim C(t) = +00: 

t^T g 

II /(*) - 9(t) Hi- + II A,(t) - X g (t) Ik- + II e 2A ' (4) - e 2A »W |k- 

+ || e^'WeJCt) - e^'We^t) |k- + || e x f^e f (t) - e x ^e g (t) \\ L ~< £(t)d 

(3.22) 

and since (|3.17|l implies the estimate: 

m,f(t,r 



r 

o2A g (t,r) 



|P/(*, 



e 2A, _ e 2A s 



W/ft, r) m g (t,r) 



+ 2 I P/ I (*' r )' 



one deduce from (|3.22|) that: 

II e/(t) - e a (t) |k- + || A/(t) - A fl (t) |k- + || /^(i) - ^(t) |k~< £(t)d 

where £ € C([0, T g [) depends only on g, strictly increasing, £(0) > 0, and 
lim £(t) = +00. Define 



e := mm < £1, 



2^(0) 
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Then a e C(]0, e[) is positive, strictly decreasing, and lirntr(/3) = T g . Let 

° o 

(/, A, e) £ D be such that < d < e. Then on the interval [0, min{er(d), T (/)}] 
the following estimate holds since £ is strictly increasing: 

Alt) := || e f (t) ~ e g (t) || £0 . + || e^e){t) - e 2A »4(t)e^(t) \\ L - 

+ II A/(i) - X g (t) + || /4(t) - /4(t) || L -< f(t)d < £(<r(d))d = -^d = i 

Thus, 

A(t) < ~, for f < o-(d). (3.23) 

Assume Ty < <r(d). Then by definition of Tb(/) and l|3.23|l we obtain the identity 
To(/) = min{Ty,T g } = Ty, in particular the estimate 

I « |< U g (a(d)) < oo, 

holds for all (ar,u) € supp/(i) and t £ [0,T/[. Since Tj < a(d) < oo, this is a 
contradiction to theorem 3.1, and we have shown that Tf > a(d). Furthermore, 
H3.23fl . implies that Tb(/) > <r(d) so that the estimates which were established 
on the interval [0, Tb(/)[ hold on [0, c(d)], and the proof is complete. 
Next, we use the above theorem and the following formula 

to obtain as it is proved in detail in [8], the 

° 

Corollary 3.1 Let (/, A, e) € D be such that 

d=\\ f-9 |U«< e. 
Then the following inequalities hold on [0, c(d)]: 

supp/(t) C {(£,«) eK 6 |i|<r +t, |«|<C^(t)} 

|| A^(i) - A^(t) + || (t f (t) - A fl (t) ||l»< S(t)d (3.24) 

i(| *,(*)-**(*) | + | A^)-A' ff (i) | + | ^{t)-^' g {t) \){t,r) < S(t)d (3.25) 

\Zf-Z g \(t,0,z)<S{t)d, for z e supp) U supp^ (3.26) 
\\T^ f (t)-T^ g (t)\\ L ^<S(t)d (3.27) 

|| R a % f (t) - R a % g {t) || £ -< S(t)d, (3.28) 

where a and e are introduced as in Theorem 3.2 and S € C([0,T ff [) is a positive 
increasing function. 
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Using once again Theorem 3.2 and Lemma 3.1 in |7|, one deduces the 

Corollary 3.2 With the assumption in Theorem 3.2, let g £ C 2 (R 6 ). Then the 
following 

II d z f(t)-d z g{t) || £ -<5(t)(||/-fl|U-+|| dj - d 2 °g \\ L ~), te[0,a(d)} 

° o ° ° o 

holds for initial data (A, e, /) £ D , with d < e, and \\ d z f — d z g \\l°^< 1, where 
e, a and d are as in Theorem 3.2 and the nonnegative, increasing function 
S G C([0, T g [) being as in the Corollary 3.1. 

4 Global existence for small initial data 

This section is concerned with the proof of a global existence theorem of solu- 
tions when the initial data and the charge of particle are sufficiently small. We 
also observe that in the corresponding spacetime, all geodesies are complete and 
then the obtained spacetime is complete. Here is the main result of this section: 

Theorem 4.1 For all ro > 0, uo > and A > there exists e > such 
that if (A, jU, e, /) is the maximal solution of the asymptotically fiat, spherically 

o ° 

symmetric EVM system with data (A, e, /), satisfying 

o o o 

supp/ c B(r ) x B(u ), II / || L oo< e, \\ A || L oo< A, 

(A, e) being a regular solution of the constraints, then the solution exists globally 
in t. Moreover, it satisfies condition (FS) stated below on R with 8—1 and 
some constant 7 > 0, 

II Pit) lUco, || p(t) || L „, II k(t) || £ «, II M(t) || £ . = 0((1+ I t I)- 3 ) 

|| A(t) |U~,|U(t) = 0((1+ [ 1 1)" 1 ) 

II r^(*) IU-, || 2V(t) IU-, || e(t) iu- = o((i+ 1 1 1)- 2 ) 

II Ks,(t) ||L^=0((l+|t|)- 3 ) 

for teM, and the geodesies are defined on R. 

Before giving a sketch of the proof of this result (reader can refer to jS] to get 
more details) , we need the following result obtained by distinguishing the cases 
r < ?'o and r > r$: 

o 

Proposition 4.1 Let (/, A, e) € D. Then the following inequalities hold: 

|e(r)|<C||/|| L - 
|£(r)|,|m(0,r)|<C||/|| L - (1+ || / |U~) 

for all r > 0, where C is a constant with depends on ro and A. and m(0, r) is 
deduced from \3.4\j setting r = 0. 
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4.1 The free-streaming (FS) condition 

Take a regular solution (/, A,/i, e) which exists on a time interval [0, T[. For 
S £]0, 1] and 7 > we consider the following decay condition on an interval 
[0,T'[C [0,T[: 



(FS) 



I A(i,r) I + I (i(t,r) I + I A'(t,r) | + | /x'(t,r) |< 7(1 + *) 
I M(t, r) I + ^ + I (t, r) |< 7(1 + t)-*- s 



-is 



i(|A(i,r)| + | M '(t,r)| + |A'(t,r)|)< 7 (l + t)- 2 - d ' 
LI e(*,r) l^lil+t)- 1 - 6 

for t e [0, T'[ and r > 0; where 

H = e- 2X + (// - A')(m' + l?j - e- 2 "(\ + A(A - £)). 

The following result, obtained by using the Gronwall inequality shows that 
under such an assumption, the momenta cannot grow very much: 

Lemma 4.1 Let 5 e]0, 1] and uq > 0. Ifj is sufficiently small and if (/, A, p,, e) 
is a solution which satisfies (FS) on an interval [0, T'[, then every solution of 
the characteristic system satisfies the estimate 

I V(t,0,x,v) \< u + 1, (x,v) e M 3 x B(u ), t G [0,T'[. 

Next, using Lemma 3.1 in [xj, one obtains as in the uncharged case, that any 
solution which satisfies (FS) allows the matter quantity p to behave like t~ 3 . To 
obtain the latter, it suffices to show that det(d v X(0, t, x, decays like t -1 

and one deduces: 



Since 



V 7 ! + v 2 f(t, x, v)dv < Cit~ 3 , t €]0, T'[. 



ie 2 ^V(i,r)<^r 



we obtain the desired result that is: 

Lemma 4.2 Let 5 s]0, 1] and ro, uq, Cq > 0. Then there exists constants 7 > 
and C\ > such that if solution (/, A, /i, e) satisfies the free- streaming condition 
(FS) on an interval [0, T'[, supp/(0) C B(ro) x B(uq) and 

o o 

II / + II ^ IU°° + II ^ + II e ||l°°< Co, 

then 

II ||i-< Ci*- 3 , te]0,r'[. 

Remark 4.1 M^e are not surprised by the fact that p decays like t~ 3 , since even 
in the Minkowski spacetime, free particles starting in a compact set spread out 
linearly with time and the associated density decays like t~ 3 as t — > 00 (see [jffi). 
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4.2 Decay of the fields 

Let (/, A, fi, e) be a solution so that (FS) is satisfied. From the following esti- 
mates 

I e(t, r)\<c\ f s 2 | M(t, s) \ ds < (7(1 + t)~ 2 
r Jo 

I e(t, r) | < / s 2 | M(t, s) \ ds < C(l + 1)~ 2 

r Jo 

I H(t,r) | < C(l+t)- 2 

one deduces that the decay of p implies the decay assumptions stated in (FS): 

Lemma 4.3 Let ro,Co,Ci > 1. Then there exists a constant 7 > such that 
if (/, A, /U, e) is a solution on [0,T[, satisfying the estimates 

II IU-< Cifi + 1)- 3 , te[o,r[, 

and 

sup{| x\\(x,v)€ su PP /(0)} < r 0) || ||l~< C 

for some T' g]0,T[, then (/, A, /j,, e) satisfies the free- streaming condition (FS) 
on the interval [0,T'[ i/ie parameters 6=1 and 7. 

4.3 Proof of Theorem 4.1 

o 

Let r , uo > and A > be fixed. Take (/, A, e) e _D, where I? is the set of initial 

o 

data introduced in Sect. 3. Using Proposition 4.1, we observe that if / is small 
in the i°°-norm, then so are ji and e. We choose e > small enough in such 

o 

a way that for all nonnegative, spherically symmetric initial data / € C^°(R 6 ) 
with 

o o 

supp/ C B(r ) x B(uo), and || / ||l=°< £ the estimates 

8ir / / + v 2 f(y, v)dvdy < r, r > 0, 
J\y\<r Jr 3 

II A ||l-< A, II / ||ioo + || fi \\ L oo + || e \\ L °°< 1 

hold. Using Theorem 3.1, we have for such initial data a local solution on some 
right maximal existence interval [0, T[ and we can choose Co = 1 when applying 
Lemma 4.2 and Lemma 4.3. Take g = X g = p g = e g = 0, and T g = 1. Applying 
Theorem 3.2, there exists e > 0, a positive decreasing function £ € C([0, 1]) and 
a positive decreasing function a e C(]0, e]) such that lirna(/3) = 1, and for any 

o 

solution (/, Xf,fif,ef) with =| / |l~< e, and the estimates below 

|| f(t) \\ L oo, II e 2A /W || L oo, || e f (t) \\ L oo< £(t)e, t e [0,a(e)] 

and then 

II ||l»< C£(i) £ , te[0,<7(e)], 
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where C = C(uq) is a constant. So, we can choose e small enough to have 

CLe < 1, where L := sup and obtain a solution (/, A, p,, e) which is defined 

te[o,i] 

on the interval [0, 1], with 

II lk»< l, te[o,i]. 

Take 5 = ^ and choose a corresponding 7 > such that Lemma 4.2 and 
Lemma 4.3 hold. Let C 7 be the constant corresponding to 7 and define 

C* := 8(C 7 + 1) 
and we take for instance T\ = ^g* + 1 to have 

7c .(l+t)- 1 < fort>Ti. 

Using Theorem 3.2 and Corollary 3.1 with g — 0, we can choose £ such that the 
solution (/, A, p, e) exists on [0, T\] and on this interval the condition (FS) with 

o 

parameters 8 — \ and 7 considered above, provided || / ||loc< e. Set 

T 2 := sup{i e [0, T[| (/, A, /x, e) satisfies (FS) on [0, £]}. 

Then by definition T 2 > T\, and using Lemma 4.2 

II p(t) \\l°o< C 7 f" 3 , te]0,T 2 [, 

and we use the fact that || p(t) ||l~< 1 for t e [0, 1], to establish the following 
inequality: 

II P (t) \\ L ~<c*{i + t)-\ te[o,T 2 [. 

Note that the inequality above is obtained, distinguishing the cases < t < 1 
and 1 < t < T 2 . Now, using Lemma 4.3, the free-streaming condition (FS) holds 
with the parameters 6=1 and ^c*i an d with the choice of T l5 (FS) holds again 
on [Ti,T2[ with parameters ^ and 6 = |. By the construction of T 2 we obtain 
T 2 = T. We deduce from Lemma 4.1 

supp/(£) Clx B(uq + 1), ie[0,T[ 

and using Theorem 3.1, we conclude that T = 00. 

Note that the decay estimates of p(t), k(t), M(t) and iV(t) come with the 
proof. So, we just have to estimate the metric, the Christoffcl symbols and the 
Ricmann curvature tensor. From 

X(t, r) = - ^ X'(t, s)ds = - J^°° e 2X ^ (-^^- + ^sp(t, s)j ds 



and 

I A(t,r) 



A(0,r)+ / \(a,r)dt 
Jo 



<|| A(0) Hioc +C I (l + s)- 2 ds 
< C 
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we deduce the following: 

rro+t ^1+ „\ r oo , f r r +t 



I A(t,r) | < C 



r ii^i ds +c r ^ds+cii+ty 3 r S d S 

JQ 3 Jra+t 3 JO 



< C'il + t)- 1 , r > 0, t > 0, 

where M > is the A.D.M mass of the solution. The estimates for /j, are similar. 
The estimate of Christoffel symbols is based on the following: 

l- e -2A(t,r) 2m(t,r) B f r ° +t . . . n 
= -^<8tt p(t,s)ds <C(l + ty 2 , 



' 1 Jo 

and the estimate of Riemann curvature tensor is deduced from the above. Now, 
from the above decay, we deduce that 

lim A(t, r) = lim l-i(t,r) = 

t — > + OG t — > + 0O 

and since we get 

II e(t) |U-< V2 || p(t) ||l-, 
one concludes that e(t, r) — > . So, the solution (/, A, (X, e) is asymptotically 

t — >+oo 

flat in time coordinate. To end the proof of Theorem 4.1, we notice that all 
geodesies are complete. To see the latter, one proceeds exactly as in the un- 
charged case (we refer to HDJ). Thus the proof of Theorem 4.1 is complete. 

5 The Restricted Regularity Theorem 

In this section, we aim to prove that a solution may be extended whenever / 

o 

vanishes in a neighborhood of the center of symmetry. Take / > 0, spherically 

symmetric, C°°, compactly supported so that (|3.1|l is satisfied. Let (X,fi, e) be 
a solution of the constraints. Then Theorem 3.1 launches the regular solution 
(/, A, fi, e) of the initial boundary value problem on [0, T] x M 6 , for some T > 0. 



Theorem 5.1 Let (/, A, fi, e) and T be as above (T finite). Assume there exists 
e > such that 

f(t,x,v)=0 if 0<t<T and \ x \< e. 
Then (/, A, fi, e) extends to a regular solution on [0, T'[ for some T' > T. 
Proof: Define 

P(t) := sup{| v\: (x,v)£ supp/(<)} (5.1) 

By Theorem 3.1, it suffices to show that P(t) is bounded on [0, T[. To do so, 
we introduce the following notation: 

u :=| v |, w := r~ l x.v, L := r 2 u 2 — (x.v) 2 . 
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Note that the square root of L is called the angular momentum. Since L is 
constant along the characteristics, we obtain: 

L = (5.2) 

r = e^ 7 =^= (5.3) 
vi+« 

w = L p^ x -wk- p!e^ x ^1 + u 2 + qe x+t "e. (5.4) 

r 3 v 1 + u 2 

Now, the values of any solution / of the Vlasov equation are conserved along 
characteristics. So, 

0</<sup/<C = C(/, £ ,T). 



Next, multiplying the Vlasov equation by vl + v 2 and integrating in v, using 
Ij3.14jl to show that 

l/2ft(e 2 V) = -qe X+tl eN/r, d t := ^ 

and the definition of A plus the relation A' + fi' = Airre 2X (p + p), yields the 
following conservation law: 

dtp + div ( e^ x [ vfdv)=0. (5.5) 

x V JR3 / 

So, the integration of (|5.5(l in i-variable yields, since / is compactly supported: 

p{t, x)dx — j p(0, x)dx = C 
and one deduces: 

< m(t, r) < p(t, £)dx <C, r > 0. 



Also, we have < L < C on the support of /. Thus 

u 2 = w 2 + L/r 2 < w 2 + C/e 2 <w 2 + C. (5.6) 

The inequality (|5.6(l shows that we need to obtain a bound for w. To do this, 
we introduce 

Pi(t) := inf{w : 3x, v G R 3 with f(t, x, v) ^ Oandw = r~ 1 x.v} 
P s {t) := sup{w : 3x, v G R 3 with /(t, x, v) ^ andw — r~ 1 x.v}. 

Since 

measure{w : (x,v) G supp/(i)} < TrCe~ 2 {P s {t) - Pi(t)), \ x |> e (5.7) 
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one deduces that if Pi{t) and P 8 (t) are bounded, so is P(t). To do so, we 
concentrate on the characteristic equation for w. Using (|3.17l) . one obtains: 

1 =e M ~ A - r" 2 me A+ Vl + u 2 + 4Trre f " +x (wk - py / 1 + u 2 ) 



rVl + u 2 ' (5.8) 

+ qe» +x e. 



As it is shown in [U], the first term in the right hand side of 1|5.8JI is bounded 
on the support of /. The same is true for the second term: 



< r- 2 me^ +x \f\ + u 2 < e~ 2 Cy/C + w 2 < C^C + w 2 . (5.9) 

We now focus on a bound of the contribution of electric field. The integration 
of l |2. 5(1 in r- variable on [0,r] yields: 

e(t, r) = -|e- A ^ r s 2 e x(t ' s) M(t, s)ds. (5.10) 
r z Jo 

The total charge Q(t) of the system is given by: 

r+ao r- 

Q(t):=4nq sV (M) / f(t,s,v)dvds. (5.10') 

Jo Jr 3 

Since the distribution function / satisfies the Vlasov equation one obtains the 
following conservation law (see Proposition 3.1, [7]): 

^ (e x f fdv) + div ( e" f —L=fdv) = 0. 

Ot \ J M 3 J X \ J R3 y/l + V 2 J 

With the help of this, we obtain: 

Q(t) = 4Trq£ s 2 J^e A f f<h)ds 

= Airq lim f s 2 ^- ( e x [ fdv) ds 
<"->+«> Jo at V Jr 3 J 

= i lim / 4" ( gA / f dv ) d y 

r ^ + °°J\y\<r dt V Jr 3 J 

= —q lim / div ( e M / - fdv) dy 

r ^ +00 J\ y \<r y V Jr 3 Vl + v 2 J 

= -Q lim / e M / fdvdu(y) 
r ^ +QO J\v\=r Jm 3 rVl + v 2 

= -Airq lim re^*'^ N(t,r), since / du(y) = Airr 2 

r^+oo J\y\=r 



0, sinceiV(t) is compactly supported. 
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Thus Q{t) is a constant and we can write: Q(t) — Q(0) := Q- So with the above 
we can deduce a bound of e(t, r), since A > 0: 

I e(t, r) |< g/(47rr 2 ) < Ce" 2 < C. (5.11) 

and we have a bound of the fourth term in the right hand side of Q5.8fl . since 
A + /i < 0: 

0<\qe^ +x e(t,r)\<C. (5.12) 
Using the above inequalities, we obtain the following estimate for w: 



-C\/C + w 2 + 47rre^ +A (wk -p^/l + u 2 ) <w< C + 4irre^ +x (wk - py/l + u 2 ) 

(5.13) 

Since < Anre^ +X < C on the support of /, we concentrate on the quantity 
wk — py/l + u 2 . Set 

w := r^^-x.v, u :=| V |, Z := r 2 {t 2 — (x.{>) 2 . 

Then 



wk — p\J 1 + u 2 = / /(i, v) ( ww — \/ 1 + ii 2 ) 
Jm? V V 1 + u 2 J 



(5.14) 



+ M 2 e 2 V. 



2 

Using once again (|5.10|l . and the fact that the total charge is conserved with 
time t one obtains: 

< e 2A e 2 < C/r 4 < Ce~ 4 < C. 
So, the above and l|5.6[) . yield: 

< 1/2^1 + u 2 e 2 V < Cy/C + w 2 < CyJC + P 2 (t). 

Next, for w > 0, if w < then the first term in the right hand side of (|5.14(l is 
negative and we obtain the following estimate for w: 



w <C + Cy/C + P 2 (t). (5.15) 

For w > 0, one has (see 

/ / (ww - xJl + u 2 ^L=) di < CP s {t) ln(l + P 2 (t)) 

JR3 V V 1 + U 2 J 

and using i|5.14[l we deduce: 



w 



<C + CP s {t) + CP s (t) ln(l + P 2 (t)) 



Set <p(x) — ln(l + x 2 ) — 1, a = \fe — 1. Then ip(a) — 0, <p(x) > for x €]a, +oo[ 
and tp(x) < for x E [0,a[. Suppose that P s (t) > a for t e [0,T[. Then 
Ps(t) < P s (t) ln(l + P 2 (t)) and we deduce from the inequality above that: 

Wl -4- P 2 (t^ 

w<C + CP s (t) ln(l + P 2 (t)) + C 1 p s(t) sU ■ (5-16) 
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If < P s (t) < a for t e [0,T[, then P s (t) ln(l + P 2 (t)) < P s (i), and we deduce 
from the above that: 

w <C + CP s (t). (5.17) 
Let w(t) be the values of w along a characteristic and set: 

to := inf{r > : w(s) > Ofors G]r,tQ. 

Then w(to) < C and (|5.16(l yields by integration: 

, ln(l + P 2 (r)) 



w(t) ^C + Cj^ [P s (t) ln(l + P 2 (r)) + v p ^ " ) dr. 

Set P s (t) := max(0, P s (t)), then for to — 0, we can write: 

Ps(t) <C + c£ \P s (r) ln(l + P s 2 (r)) + Ml t ^ {t)) ^ rfr 

and using the generalized form of the Gronwall inequality (see |15j . page 26), 
one has: 

P.(t) <z(t), 

where z is a solution of the following differential equation: 



z(t) =Cz{t) ln(l + z 2 (i)) + C- 



z(t) 



and since z(t) = \J exp(e 2C *) — 1 is a solution of the equation above, one finally 
obtains: 



Ps{t) < P s {t) < v /cxp(e 2 «) - 1 < C. (5.18) 
Also, we integrate (I5.17f) on [0, t] and obtain: 

ft 



P s {t)<C + C j P s {r)dT 
Jo 



and deduce by Gronwall's lemma the following bound of P s (t): 

P s (t) < max(a,e c *), te[0,T[. 

Note that the above is valid in the case P s (i) < 0. We now look for a bound 
of Pi(t). Suppose Pi(t) < and consider w in supp/ with Pi(t) < w < 0. For 
w < 0, we get (see [Hp: 



WW 



Vl + u 2 ^= > -C- 



y/1 + V? \A + W 2 

and one uses l|5.18|) to obtain, for w < < w < P s (t): 



ww - yl + u 2 ^ >Cw- cVc + w 2 . 

VI + 71 2 
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So, we deduce: 

/ / (ww - y/l + u 2 ^£=) dv > -CP?(t)—!-—+CP a (t){w-y/C + w 2 ). 
Jr3 V vl +vrj Vl + w 

w > -Cy/C + w 2 - CP? (t) -= + Cw 

Vl + w 2 

w 2 = 2ww <C + CP 2 {t). 

Set 

ti := inf{r > : w(s) < for s£]t,([}. 
Then > > — C and one proceeds as before to obtain: 

P?(t)<C + C f 'Pi(r)dr, P i (t)<0 ) 
Jo 

and the Gronwall lemma yields: 

F?(t) < Ce ct < C, fori e [Q,T[. 
Besides, if P»(t) > then 

0<P<(i) <P,(t) <C, 

and Pj(i) is also bounded for that case. So, P(t) is bounded and Theorem 5.1 
is proved. 

We now focus on what happens when we are far from the center. 



6 The Exterior Problem 

Let r\ and T be positive real numbers. Consider the exterior region: 

W(T,n) := {(t,r) : < t< T, r>n+t}. 
In what follows, the initial value problem (|2~3|) . JO), (f2~5j) . lf2~§jl . and 

o 

<|2.10[l are going to be studied on a domain of this kind. Take / > 0, spheri- 
cally symmetric, C°°, compactly supported and defined on the region | x |> r\ 

and let (A, /i, e) be a regular solution of the constraint equations on that re- 
gion. Since we are away from the center of symmetry, we have to change 
the boundary conditions (|2.9[) as follows: Let > Air J r s 2 p(0, s)ds, and 

o 

Moo > J r 7 s 2 e x ^M(0,s)ds. Take any solution of (JO), and lJ53)l 

on W^T, ri) and define 

fm(t,r) := rrioo - 4^^°° s 2 p(t, s)ds 

1 M(t, r) := - J, 00 S 2 e A (*' s )M(t, s)ds ( ' j 
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Then for all r > ri, m(0,r) > 0. So, we replace (12.91) by: 

e -2A(t,r) = 1 _ 2 m(ti). lim X (t,r) = lim /*(*, r) = 
e(i,r) = 4e~ A(t ' r) M(i,r); lime(f,r)=0 

Note that the restriction of a solution of the original problem with boundary 
conditions (|2.9|l on W{T,r{), satisfies (|6.2[) if is chosen to be the A.D.M 
mass of the solution on the entire space, and M M = Q/(4irq), Q being the total 
charge given by (5.10'). Besides, to obtain a local existence theorem as in [7], 
we use the fact that (I6.2[l must hold on the initial hypersurface t = that leads 
to the following condition: 

(/ \J\ + v 2 °f(x,v)dv ] dx < r -, r>n (6.3) 

c|>r 



Theorem 6.1 (Local existence and uniqueness) Fix moo > 0. Let 

o 

f > be a spherically symmetric function on the region \ x \> r± which is C°° 

and compactly supported. Let (A, fi, e) be a regular solution of the constraint 
equations on the region r > r± ■ Suppose that holds and that 

(^J + v 2 f(£, v)dv^j dx < moo (6.4) 
e X{i) (J °f(x,v)dv^j dx < AttMoo- (6.4') 

|5|>ri 

Then there exists a unique regular solution of \2.'A) . ]2.°J\) . and \2.5}) on the 

region W(T, r{) with (/, A,/i,e) which satisfies the boundary conditions ^6.2}) . 

o 

Proof: Let / be as in Theorem 6.1. Suppose 

supp) C B(R) x B(R'), R, R' > 

where B(R) is the open ball of M 3 with radius R. We only have to check that 
one can construct a solution satisfying the constraints and use to obtain the 
desired result. If r\ > R, then the constraint equations reduces to the static 
Einstein-Maxwell system and a solution for this system on the region r > r\ is a 
part of the Reissner-Nordstrom solution 0]. Now, if r\ < R, since the constraint 
equations is invariant with translation we can use perturbation techniques to 

o 

establish a solution for the corresponding Cauchy problem provided / satisfies 
an appropriate condition like l|3.1() . The reader can refer to [Hj, to have more 
details. So Theorem 6.1 is proved. 
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We now derive a corresponding continuation criterion, with proof being as 
that of Theorem 3.2 in Before doing so, we mean by maximal interval 

of existence for the exterior problem, the largest region W(T,r\) on which a 
solution exists with given initial data and the parameters r\ and moo. 

o 

Theorem 6.2 (Continuation Criterion) Fix m,oo > 0. Let f > be a 

spherically symmetric function on the region | x |> ri which is C°° , compactly 

supported and satisfies \6.S\) and \6.4\j - Let (A, fx, e) be a regular solution of the 
constraints on the region r > r-± so that (6.4') holds. Let (f,X,fj,,e) be a regular 

solution of if 2.21) . if£.3j) . \2.$ and if 2.5)) on W{T, r\) with initial data (/, A, ft, e). 

If T < oo and W(T, r\) is the maximal interval of existence, then P is bounded. 

7 The Regularity Theorem 

Here we consider the initial boundary value problem l|2.2fl . (|2.3f) . (|2.4fl . 12. 5f) . 
I|2.9fl and H2.10f) . The following result shows that if a solution of this Cauchy 
problem develops a singularity at all, the first one must be at the center. 

Theorem 7.1 Let (A, fi, e) be a solution of the constraint equations. Let (/, A, [i, e) 
be the corresponding regular solution defined on a time interval [0, T[. Suppose 
that there exists an open neighborhood U of the point (T, 0) such that 

sup{| v |: (t,x,v) e supp/n (U x K 3 )} < oo. (7.1) 

Then (/, A, /i, e) extends to a regular solution on [0, T'[ for some T' > T. 

Proof: See [TT| . 
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